Dephasing representation of quantum fidelity for general pure and mixed states 
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General semiclassical expression for quantum fidelity (Loschmidt echo) of arbitrary pure and 
mixed states is derived. It expresses fidelity as an interference sum of dephasing trajectories weighed 
by the Wigner function of the initial state, and does not require that the initial state be localized in 
position or momentum. This general dephasing representation is special in that, counterintuitively, 
all of fidelity decay is due to dephasing and none due to the decay of classical overlaps. Surpris- 
ing accuracy of the approximation is justified by invoking the shadowing theorem: twice-both for 
physical perturbations and for numerical errors. It is shown how the general expression reduces to 
the special forms for position and momentum states and for wave packets localized in position or 
momentum. The superiority of the general over the specialized forms is explained and supported 
by numerical tests for wave packets, non-local pure states, and for simple and random mixed states. 
The tests are done in non-universal regimes in mixed phase space where detailed features of fidelity 
are important. Although semiclassically motivated, present approach is valid for abstract systems 
with a finite Hilbert basis provided that the discrete Wigner transform is used. This makes the 
method applicable, via a phase space approach, e. g., to problems of quantum computation. 
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I. INTRODUCTION 

Time evolution in classical mechanics is very sensitive 
to perturbations of both initial conditions of a trajectory 
and the Hamiltonian. Because of the unitarity of quan- 
tum evolution, on the other hand, the overlap of two 
different quantum states remains constant in time. How- 
ever, we can still define sensitivity of quantum evolution 
to perturbations of the Hamiltonian. This is usually done 
using the notion of quantum fidelity (sometimes called 
Loschmidt echo), defined for pure states as pj 



M(t) 



e +iHH/n e - 



iH°t/h 



(1) 





Here \t/j) is the initial state, H° and H e — H°+eV are the 
unperturbed and perturbed Hamiltonians, respectively. 
In words, fidelity is the overlap at time t of two identical 
initial states evolved with two slightly different Hamilto- 
nians. Because of its relevance in theories of decoherence 
and in experimental realizations of qua ntuiri ^nuo t rta- 

as well as classical 
fidelity has been extensively stud- 
ied in the last few years. Many universal regimes of fi- 
delity decay have been found in different limiting cases 
SSSSSSSIIimEi. Many of these works 
used a semiclassical approach, but before Ref. only 
as a starting point for various approximations, because 
of difficulties in treating an exponentially growing num- 
ber of terms in the general semiclassical expression for 
fidelity, especially in chaotic systems. This problem was 
solved in Ref. [13| by a uniform expression for fidelity 
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which implicitly summed over all these contributions us- 
ing an integral over initial conditions, similar in spirit to 
Miller's initial value representation |43.l4^ . This surpris- 
ingly simple and accurate expression, although limited to 
wave packets localized in position, has been successfully 
applied as a starting point to derive fidelity decay in the 
deep Lyapunov regime 1351 and the plateau of fidelity 
in neutron scattering |l2j . Five other known regimes of 
fidelity can also be simply described by this method 15]. 

In a recent Rapid Communication |l4j , the uniform ex- 
pression for fidelity was justifie d by the shadowing the- 
orem of classical mechanics 0, 0| and a more general 
and, in fact, always more accurate expression, valid for 
arbitrary pure states, was stated. One purpose of the 
present article is to provide (in Sec. [Hj a detailed deriva- 
tion of this general semiclassical expression (|18H for fi- 
delity of arbitrary pure, i. e., also nonlocal states. Fi- 
delity is expressed as an interference sum of dephasing 
trajectories weighed by the Wigner function of the ini- 
tial state. The general derivation provides an alternative 
and more explicit justification of the validity of this de- 
phasing representation (DR). Interestingly, in Sec. IIIII 
it is shown that the same dephasing representation is 
valid also for general mixed states. Section IIVI shows 
how the general expression reduces to the original form 
|13| and other specialized forms for position and momen- 
tum states or Gaussian wave packets localized in position 
or momentum. In Sec. [V] the general dephasing repre- 
sentation is tested on a non-local state-a coherent super- 
position of two separated wave packets-and on two two 
types of mixed state-an incoherent superposition of two 
wave packets and a completely random state. It is also 
shown that the general expression is superior to the orig- 
inal form [l3| even for a single Gaussian wave packet. 
All numerical calculations are done for a system with 
a finite Hilbert basis. In such systems, quantum phase 
space can be rigorously defined if the original Wigner 
function (4^ is replaced by the discrete Wigner trans- 
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form [43, |4!J, |5(J, |5lj, |52( • Since this discrete transform 
can be defined in a general abstract Hilbert space with 
finite basis, the present approach should be applicable 
to problems of quantum computation if phase space ap- 
proach is used [53j. In Sec. IVII DR is compared to other 
"Wigner" methods. The main conclusions of the paper 
are summarized in Sec. IVIII 



II. DEPHASING REPRESENTATION FOR A 
GENERAL PURE STATE 

Fidelity amplitude for a general pure state \tp) can be 
written as 



0(t) = ^\e +lHH / h e- m "^ h U) 



(2) 



In order to derive the general dephasing representation 
of fidelity, we could start by replacing the two quantum 
propagators in Eq. Q by the corresponding semiclassical 
Van Vleck propagators |H3|, as in Refs. 0,13 • However, 
we will save some effort if we start directly from the semi- 
classical initial value representation (IVR) |43ll44| for the 
two Van Vleck propagators, 



-iH°t/h 



-d/2 



dr' / dp' \dr' t (r' ,p' )/dp' 



,1/2 



A. Uniform semiclassical expression for fidelity 

If we further expand the 5 function in integral © as 
an integral over a dummy momentum q, 

(rJV t > = 6 (Ar t ) = (2nhy d J rfqe^^^, 

we obtain a "full" uniform semiclassical expression for 
fidelity, 



xe ^(r^; t )M |r , )(r , L (3) 
e+ WH/H „ ( _ 27rffi) --/2 j K J dp « \ dr >> p //) /dp >>f/* 



x e -^ e (ro,Po ;0/ ft | r ^}(r"l 



Here rg,pg and r ',p ' are the initial conditions of tra- 
jectories of H° and of H e , respectively, and r' t ,p' t and 
r" , p" are the corresponding coordinates and momenta 
at time t. Action S ' of a trajectory of the unperturbed 
Hamiltonian H , is given by 



Ou„if(t) = {2nh) 



-2d 



dr' J dp' J dr'ij dp'^J dq 



dr' t 


1/2 


dr' t ' 






dPo 



1/2 



r (r ')^(r ) 



x exp \ - [S ' - S e " + q • Ar t ] 



(0) 



This integral is, formally, semiclassically "exact." In par- 
ticular, it is free of caustics, unlike, e. g., the Van Vleck 
propagator. Because it is expressed only in terms of ini- 
tial conditions (and dummy momentum q) , it appears to 
be ready for numerical evaluations. Unfortunately, this 
integral is highly oscillatory, and very difficult to com- 
pute, especially in many-dimensional or chaotic systems. 
Therefore we will take an alternative route, using a fur- 
ther approximation, but obtain an integral much easier 
to tackle numerically. 



B. Dephasing representation 

First, let us make a change of variables 
{r',r",p',p"} — > {r, Ar,p, Ap} in integral 10. It 
should be emphasized that we do not assume Ar or Ap 
to be small. New variables (averages and differences) 
are defined for all times from to t as 



r = ±(r' + r") 



(7) 



S°' (r' 0) p' ;t) = f dr [p' T -r' T - H° (r' T , p' T ; r)] 
Jo 



(4) 



P=2(p' + p"), 
Ap = P " P ', 



Similar expression holds for the action S e " (r ',p ';i) of 
a trajectory of the perturbed Hamiltonian H t . In the 
simplified notation above, the square roots of the deter- 
minants in Eq. (J3J) also include the appropriate Maslov 
indices [55|. Using the IVR expressions ©, fidelity am- 
plitude J5J) becomes 



Oivr (t) = {2irhy d / dr' / dp' 



dr' ' 



dp'o' 



dr' 



dp'o 



1/2 



dr" 



dp'o' 



1/2 



<V>|r ')<r;'|r;><r Q |V>> 



J(S°'-S*")/H 



(5) 



The Jacobian of this transformation is unity. If we in- 
tend to perform integrals over Ar and Ap first, we can 
consider rg and po as fixed for the moment, and write 



dr' t 


1/2 


dr'l 


1/2 


d (-Ar t ) 


1/2 


<9Ar t 


1/2 


dAr t 


dp'o 




dp'o' 




d(-Apo) 




dA Po 




cApo 



(8) 



0(t) = (2tt/j) 



dr / dp / dAr / dAp 



dAr t 



x rK)S(Ar t )^(r' ) exp 



dAp Q 

(9) 
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Next we change variables from Apo to Ar t and eliminate 
the 5 function, 



0(t) = (2nhy d / dr / dp / dAro^* (r ') tf> {r' ) 



x exp 



5 £ ") 



(10) 



Ar t =0 



The present form is equivalent to Eq. @. On one hand, 
the present form appears much simpler (a 3d- vs. 5c?- 
dimensional integral), on the other hand it is not an in- 
tegral over independent variables because it contains a 
constraint on the final positions (Ar t = 0). 

While we do not intend to evaluate this integral by the 
stationary phase (SP) approximation, it is instructive to 
check where the action difference S ' — S e " is stationary 
because those regions give the main contributions to the 
integral. Variation of action S 0/ gives 

6S Q ' = -p' -5r' +p' t -5r' t 

and a similar expression holds for SS e ". Due to the 
Art = constraint, we have a constraint Sr' t — Sr" on 
the variation of endpoints, and therefore 



5 S' 



01 



S' 



-Po 



Sr' 



Po 



K 



A Pt • 6r' t . 



Expanding variation 6r' t in terms of variations 6r' and 
Sp' , we find 



S(S° 



Apo - Apt 



dr' t 



dr' 



dr' 

Apt • 3rf 



Po 



K 



5Ar . (11) 



Note again that so far we have not assumed anything 
about closeness of the two trajectories. Since we can eas- 
ily shift integration variables ro and po to r and p' in 
Eq. I |10|l . variation (|ll|l indeed tells us where the action 
difference would be stationary. There are three station- 
ary phase conditions, 



dv' t 



A Po - Apt • = 0, 
Br' 

p ' • <5Ar = 0. 



(12) 

(13) 
(14) 



The third SP condition was intentionally written in the 
full form. In general, all three conditions would be satis- 
fied only for a discrete set of trajectories (3c? equations for 
3c? unknowns). However, if the perturbation were e = 0, 
one could immediately guess that there is one continu- 
ous set of solutions satisfying Apo = Ap t = Aro. The 
first two conditions are satisfied exactly, the third one 
approximately for small variations ^Ar . Even though 
the third condition is satisfied only approximately, we 
obtain the correct result -identical trajectories Ar r = 



for all times t, < r < i-and as we shall see below, 
also the final result for fidelity will become exact in this 
limit (e = 0). If we add the perturbation, these precise 
solutions break down, due to the exponential sensitiv- 
ity of classical dynamics. However, as was shown in Ref. 
[Tj|, if the shadowing theorem |43,|4£j is applicable in the 
given system (for a given perturbation e and up to time 
t) , there will be a very near solution with Ar T w for all 
times r, < r < t. Putting off a discussion of the shad- 
owing theorem until later, suffice it to say that this the- 
orem, completely counterintuitively, guarantees that we 
can compensate one exponential sensitivity (to perturba- 
tions of H ) by another exponential sensitivity (to initial 
conditions) and get a trajectory which remains very close 
to the unperturbed trajectory up to time t. In fact, these 
approximate ( "diagonal" ) solutions with Ar T w will be 
by far the most dominant ones because for short times 
no other solutions exist and for long times the diagonal 
solutions dephase much slower than the remaining ( "off- 
diagonal") solutions with different trajectories. Again 
this will be justified later in this section. Assuming the 
validity of shadowing, the "diagonal" solutions dephase 
as 

5°' - S'" we / drV (r T , r) - Ar t • p t + Ar • Po (15) 
Jo 

= -ASj - Ar t • p t + Ar • p . (16) 

The first term is due to the perturbing potential eV along 
the unperturbed trajectory, the other two terms are due 
to the small difference of trajectories at time t and at 
time 0. Substituting this action difference into integral 
Hltty . we obtain the dephasing representation 

Pdr(*) = (2nhy d J dr J dp J c?Ar 

x ijj* (r + iAro) i/) (r - ^Ar 



- (-AS t + Ar • po) 



x exp 

The final result is more succinctly written as 

Odr (t)= J dr Q J dp p w (r ,p ) exp(-iAS t /K) 
using the Wigner function of the initial state \ijj), 
pw(r,p) = (2irfrr d J dAryj* (r + ^Ar 



(17) 



(18) 



x tp ^r — 2^ r J exp ' ' 

The general expression Qlfcty expresses fidelity as an 
interference integral over initial positions ro and mo- 
menta po- Because of this property it was called dephas- 
ing representation in Ref.jlj]. The amplitude of each 
term is given by the Wigner function pvv (ro, Po) and the 
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phase by the integral of the perturbing potential along 
the unperturbed trajectory, ASt ( r OiPo)- This is a very 
intuitive and simple picture that differs from the sim- 
plest "semiclassical" picture only in using the Wigner 
function instead of the classical phase space distribution 
Pclass (ro, po). 

For zero perturbations, e = 0, expression l|18|) correctly 
reduces to the obvious exact result, 

Odr (*) = / dt J dp Q p w (r , Po ) = 1 (20) 

for all times t, where the basic property of the Wigner 
function was used. 

Although we started our derivation for a pure state, 
we ended up with a dephasing representation in terms 
of the Wigner function. Since this function can also be 
defined for mixed states, it appears that expression l|18|) 
should remain valid for mixed states, with appropriate 
generalization of the notion of fidelity. In Section ITTT1 it 
will be shown that this is indeed the case. 



cases in Sec. IIVI Wigner function is particularly sim- 
ple for position and momentum eigenstates (just a delta 
function), for Gaussian wave packets (a Gaussian in both 
position and momentum), or for a random mixed state 
(a constant over the whole phase space). These distribu- 
tions can be easily sampled using standard methods. For 
general pure or mixed states one can resort to a Monte- 
Carlo procedure, e. g., using the Metropolis algorithm, 
which is frequently done for the IVR approximation |44| . 

One might object that numerical computation of tra- 
jectories, due to the exponential sensitivity of classical 
evolution, will destroy the validity of the DR (|18|l . How- 
ever, here the shadowing theorem helps again-in fact in 
its original form |4(| where the perturbation was in- 
deed due to errors of numerical p ro p agation. The shad- 
owing idea, as stated in Refs. |45l l46j guarantees that for 
each numerical (noisy) trajectory there will be a nearby 
exact trajectory of H°. 

D. Comparison of diagonal and off-diagonal terms 



C. Shadowing theorem and its double use 

1. Trajectories of H° and 

Shadowing theorems in general state that (under cer- 
tain detailed conditions) for small enough e there is a time 
t such that for a trajectory of H° with initial condition 
r' , p' Q there exists a trajectory of H e with initial con- 
dition r ', Po remaining within a certain small distance 
from the first trajectory up to time t. In uniformly hyper- 
bolic systems this shadowing time t is infinite |56l l57j , in 
more general systems at least finite |4^. Since it is very 
difficult to find the maximum shadowing time t and the 
corresponding bound on the closeness of trajectories for 
a specific system, the derivation of dephasing represen- 
tation of fidelity assumed that shadowing was applicable 
for a given perturbation and time: the numerical results 
will provide the final verification. 



2. Numerical evaluation 

In order to use DR in numerical applications, one only 
needs to generate initial conditions ro, po from a distri- 
bution given by the Wigner function pw, run trajectories 
with the unperturbed Hamiltonian H° and compute the 
action difference ASt = — e / drV (iv,t) along this tra- 
jectory. There is no need to compute Van Vleck deter- 
minants or Maslov indices as in many other semiclassical 
applications. Because the Wigner function, unlike classi- 
cal probability, can be negative, some care must be taken 
to sample from its distribution. The simplest possible 
recipe would be to sample according to the probability 
\pw\ and attach a sign afterward together with the de- 
phasing factor. As we will see from the analysis of special 



Let us attempt to quantify the validity of the DR by 
comparing the importance of diagonal and off-diagonal 
terms in fidelity amplitude. These should be distin- 
guished from the "diagonal" and "off-diagonal" terms 
in the fidelity itself (i.e., the amplitude squared), which 
have been frequently discussed in the literature, see, e. 
g. Refs. 0,0 where the off-diagonal terms in the fidelity 
amplitude are already neglected. 

For short enough times t, it is clear why DR is 
accurate-there will be no off-diagonal contributions be- 
cause there will be no off-diagonal SP solutions of Eqs. 
(11211 . For long times, the number of off-diagonal solutions 
increases, but in the semiclassical limit (small K) and for 
small perturbations e, their contribution is again negli- 
gible, due to their much faster dephasing. Let us see in 
detail how this happens. 

If the unperturbed potential is denoted W, then the 
off-diagonal solutions dephase as 

ASoff-diag = S°' - S°" = f dr[W(r' T )-W(r';)} (21) 

Jo 

because for small enough e, the perturbation V is really 
unimportant in dephasing of off-diagonal terms. Assum- 
ing for simplicity that the off-diagonal terms have the 
same weight and that their action differences are Gaus- 
sian distributed, their average will be 



cxp -((AS* off _ diag ) 2 )/2fi : 



(22) 



In chaotic systems, ^(AS^ff-diag) J = 2Kyyt, where the 

diffusion coefficient is K\v = J °° dt Cw (t) and Cw is the 
potential correlator, C w (t) = (W [r (*)] W [r (0)]). 

Similar analysis can be done for the diagonal terms 0, 
HE1. Their average is then given by a formula analogous 
to Eq. , except with ASdiag given by expression l(T5|l . 
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The variance is now given by < (AS, 



2Kye 2 t. where p° (t), p e (t) are the evolved density operators, 



Because the diagonal contributions are weighed by the 
Wigner function, their total contribution is roughly equal 
to the average. The number of discrete off-diagonal semi- 
classical contributions should for long times grow as e 7 * 
where 7 is the topological entropy. Then in the worst pos- 
sible scenario, where each off-diagonal term contributes 
by its full weight (as if the Wigner function-in the case of 
diagonal terms-were unity everywhere), the ratio of the 
sum of the off-diagonal contributions to the total contri- 
bution of the diagonal terms should be 



off-diag. 
diag. 



exp { [- (K w - K v e 2 ) /h 2 + 7] t} 



For small enough e and small enough h, the off-diagonal 
terms will become negligible. Namely, the diagonal terms 
will give a smaller contribution if both h 2 < K\y /j and 
e 2 < (K w - h 2 -f) /K v . 

Similar analysis is possible for integrable systems 
0, E|- There the number of off-diagonal contribu- 
tions grows only algebraically, ~ t a and variance of 

their action difference ((AS ff-diag) 2 ) = C^t 2 where 



= linit^oot 1 J dr Cw (t). Similarly, for diagonal 
terms, ^(AS diag ) 2 ^) = C^t 2 e 2 7}. In this case, the ratio 
of the two types of contributions is 

off-diag. 



diag. 



t a exp [- (C, 



00 /^ioo 
W 



t/ (2H 2 )] 



and the condition for negligibility of the off-diagonal 
terms in the limit t — > 00 is e 2 < C^/Cy. 



III. DEPHASING REPRESENTATION FOR A 
GENERAL MIXED STATE 

There are several ways to generalize the pure-state def- 
inition (J2J of fidelity to mixed states. The simplest gen- 
eralization is 



O (t) = tr (e 



iH°t/h +iH € t/h 

pa 



(23) 



where p is the density matrix of the mixed state, normal- 
ized such that trp = 1 7]. For pure states p = \ijj)(ip\, 
this general definition reduces to the pure-state defini- 
tion (J2J. One interpretation of the general expression 
(|23[) is that the ket vectors evolve with the unperturbed 
Hamiltonian H° and the bra vectors with the perturbed 
Hamiltonian H € . Another interpretation is that expres- 
sion lll'-'il) is simply an average of fidelity amplitudes of 
pure-state components of the given mixed state. This 
should be distinguished from the often studied averaged 
fidelity 

The second possible generalization of the notion of fi- 
delity to mixed states replaces the expression for fidelity 
JTJ, rather than fidelity amplitude by an expression 

M(t)=tv[p°(t)p^(t)] =tr[/>(0)p(i)], (24) 



_ e -iH c t/h pe +iHH/h 



or, alternatively, p(t) is the evolved operator 

p(t) = e +*H't/h e -iH°t/h pe +iH°t/h e -iHH/h^ 

Again for pure states p = definition (|24|l reduces 

to the pure-state definition (|TJ. 

Finally there is another, more intuitive but also more 
complicated generalization, which uses the notion of "pu- 
rity fidelity" -the trace of the squared reduced density 
matrix |21| . 



P F (t) = tr trp(i) 

S IE . 



(25) 



where subscripts E or S denote that the trace operation 
is performed on the environment or system degrees of 
freedom, respectively. For details see Ref. [2l|. Purity 
fidelity (|25|l does not, of course, reduce to the definition 
of fidelity for pure states QJ. 

While dephasing representation expressions are possi- 
ble for the last two generalizations, in what follows the 
simplest generalization (|23|l is assumed. With the mixed- 
state definition the semiclassical derivation in Eqs. 
(I5|)- (lf 8} , can be followed closely for mixed states, if we 
replace the product (V>|io')(rg|7/j) = -0* (r^) ip (r' ) in Eqs. 
©, and l(T7|l by the matrix element (rg|p|ro) of the 
density operator. For instance, Eq. 10 will become 



0(t) = (27rft)" 



dr / dp / c?Ar n / dAp 



dAv t 



dAp 



xS(Ar t ) (r' \p\4) exp 



At the end, we obtain the same final result (|f 8fl . only the 
Wigner function of a pure state l|f 9|l must be replaced by 
the Wigner- Weyl transform of the density operator, 

pw (r,p) = {2TThy d JdAr + ^Ar |p| r - i Ar^ 
x exp(iAr ■ p/K) . (26) 

IV. SPECIAL CASES 

For a position state |R), ip (r) = S (r — R), the Wigner 
function l|f 9|) is 



/4 OS (r, P) = (2nhy a / dx 5 r + -x - R 



x S — -x — Rj exp (ix • p/K) 
= (2Trhy d J(r-R). (27) 
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Substituting Eq. i|19|) m to the general dephasing repre- 
sentation (|18fl . we find 

0&°R St ' ( 4 ) = ( 2 ^)- d /" d Po exp [-iASt (R, po) /ft] , 

in agreement with Eq. (1) from Ref. and with 0, 

For a momentum state |P), "0 ( r ) = 
(2nh)~ d ^ 2 exp (iP • r/ft), the Wigner function l(T§|) 
becomes 

Pw° m - St ' (r, P) = (27rft)- 2d /" dxexp [i (p - P) • x/ft] 

= (27rft)- d J(p-P), (28) 
and the general DR of fidelity (|18|l reduces to 



For general (non-Gaussian) wave packets, which are nev- 
ertheless localized cither in position (about R) or mo- 
mentum (about P), we can use the general property of 
the Wigner function 



drp w (r,p) = |V>(p)| , 
y dpp w (r,p) = \ip(r)\ 2 , 
and obtain, upon substitution into the general DR <|18l) . 



PO rS .w. P . (t) = J (R) /h] ^ (po)| 2 



, 0~ p ' (t) - / dr exp [-iA5 t (r 0) P) /ft] \i> (r„) 

O mom.st. (f j = {2nh) -d J ^ (r0) p) J 



(32) 



(33) 



A general Gaussian wave packet with average position 
R, average momentum P, and position spread a, 



-d/4 



tp (r) = (7rcr 2 ) ' exp 
has Wigner function 



iP-(r-R) /ft- (r-R) 2 /2<r 2 



p^-G-w.p. (r . p) 



[ 7T(7 ) (27m) 



c/x 



x exp <; - (p 



(r-R) 2 + (x/2) 2 /a 



(irhy a exp -(r-R)V^-(p-P) cr 2 /ft 

'(29) 

In general the dephasing representation of a Gaussian 
wave packet is (|18|) with the Wigner function l|29ll where 
we must include dephasing trajectories with varying both 
positions and momenta. Only in special cases, such as 
when the wave packet is strongly localized in position 
(i. e., when a <C ft 1 / 2 ), can we make a further simplifi- 
cation by replacement of AS t (r ,po) by AS t (R, p ) in 
Eq. I|18l) . Then we can compute the ro integral in Eq. 
(|18|l analytically and obtain 



O 



pos.G.w.p. 
DR 



(t) = (aVvrft 2 ) 



dp Q 



x exp 



-iASt (R, Po) /ft - (P - P)' ^ /ft J . (30) 

in agreement with Eq. (8) in Ref. 01 . There the same 
result was obtained by linearizing the Van Vleck semiclas- 
sical propagator about the central trajectory. In Section 
IVl it will be shown that the symmetric expression 129|) 
based on the general DR (|18|l is superior to the special- 
ized form (|30fl . Similarly, if the initial Gaussian wave 
packet is localized in momentum (i. e., when a 3> ft 1 / 2 ), 
we can replace ASt (ro, Po) by ASt (ro, P) and obtain 



Qmom.G.w.p. 



(*) = M 



dr 



x exp 



-iAS t (r ,P)/h-(T-R) 2 /a" 



(31) 



Finally, for a completely random state, i. e., an inco- 
herent superposition of all pure basis states, the density 
operator as well as its Wigner function l|26|) is just a 
constant (independent of position or momenta), and for 
a system with a finite phase space volume f2, the DR 
becomes 



O random st. 
DR 



(t) = ^ y dr J dp exp {-iAS t /K) . (34) 



It should be pointed out that while names like "posi- 
tion" or "momentum" states have been used to describe 
the special cases, they do not necessarily need to be eigen- 
states of the usual position or momentum operator. In 
the case of abstract Hilbert space with a finite basis, "po- 
sition" states are simply the basis states (called computa- 
tional states in the setting of quantum information, could 
be, e. g. spin eigenstates) , and "momentum" states are 
simply the states defined by the discrete Fourier trans- 
form of the original basis states |5^| . In Ref. j5^| , this 
generalized phase-space representation is used to show 
that for quite a few interesting operations on computa- 
tional states, the Wigner function evolves classically. In 
all these cases, the dephasing representation described in 
Sees. JlllVl should be applicable if discrete Wigner func- 
tion [53j is used and simple other modifications are made 
to account for the finite-size of phase space. In fact, this 
is done in the numerical examples in the following sec- 
tion. 



V. NUMERICAL TESTS 

Now let us apply the theoretical analysis from previ- 
ous sections to a specific system, the Chirikov standard 
map. Its advantage is that it is discrete, coordinate space 
is only one-dimensional, but at the same time standard 
map already contains generic complexities of classical dy- 
namics. Specifically, the phase space is mixed and so 
various simplifications applicable in quasi-integrable or 



strongly chaotic systems are in general not applicable. 
Standard map is a symplectic map defined on a compact 
two-dimensional phase space-torus, as follows, 

Qj+i = Qj + Pj ( m °d 27r) 

Pj+i = Pj - w ' (Qj+i) ~ eV ' fe'+i) ( m °d 2tt), 

where q and p are position and momentum on the torus, 
potential W (q) — — fccosq, and the perturbation is 
V (q) = — cos 2q. Using an n-dimensional Hilbert space 
for the quantized map fixes the effective Planck constant 
to be h — (2ttii)~ . (We are using letter q for the co- 
ordinate to distinguish this special system from the gen- 
eral considerations. Similarly, we will use letter Q to de- 
note the position of a position state or center of a wave 
packet.) Parameter e controls the strength of perturba- 
tion. For e <C 1, the map is close to being integrable, 
for e 3> 1, the map is strongly chaotic. The goal of this 
section is not to use the dephasing representation to ex- 
plore various universal regime that occur in these two 
limits and have been carefully studied in the literature. 
This was already done in Refs. 0, 0. The goal of 
this section is rather to explore the detailed features of 
fidelity in non-universal regimes. The optimal region of 
parameter space is in the vicinity of e = 1, since there 
phase space has a significant amount of chaotic as well 
as integrable regions. Mixed phase space is in general 
the hardest to treat and therefore this setting is chosen 
here because it provides the most challenging test for any 
approximation. 



A. Gaussian wave packets 

One might think that the general dephasing represen- 
tation (|18|l is only useful for highly non-local states and 
that the original expression 1)27(1 from Ref. E3 is good 
enough at least for Gaussian wave packets. This subsec- 
tion demonstrates that even for Gaussian wave packets, 
the general dephasing representation (|18H is superior to 
the original expression 1)2 7)1 from Ref. • 

Figure 1 compares three approximations to compute 
fidelity of Gaussian wave packets with the exact result: 
the expression 

(22) from Ref - El 

for wave packets local- 
ized in position (red dashed line), corresponding expres- 
sion (|3 ll> for wave packets localized in momentum (blue 
dotted line), and the general DR 1)18)) . with the Wigner 
function 1)29)1 . symmetrically treating position and mo- 
mentum (black solid line). The exact fidelity, computed 
by exact quantum evolution using the Fast Fourier Trans- 
form algorithm, is represented by solid dots. The param- 
eters are n — 1000, k — 0.95, e = 0.015, and the wave 
packet is localized at Q = 0.7ir and P = OAir. The num- 
ber of classical trajectories used in the calculations is 
1000. Wave packets used in parts a), b), and c) of Fig. 1 
have position spread a equal to 0.0047T, 0.167T, and 0.047T, 
respectively. For a wave packet localized in position in 
Fig. la), the original expression (|27|l from Ref. 13| works 
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FIG. 1: Fidelity for a Gaussian wave packet centered at 
Q — 0.77T, P — 0.4-7T in a perturbed standard map with 
n = 1000, k — 0.95, e = 0.015. Comparison of the exact result 
(solid dots), general dephasing representation (black line) and 
its specialized forms for position-like (red dashed line) and 
momentum- like wave packets (blue dotted line). The initial 
position spread a of the Gaussian is: a) 0.004-7T, b) 0.16-7T, and 
c) 0.047T. The DR calculations used 1000 classical trajectories. 



very well and is almost indistinguishable from the general 
DR (TSJ, as expected, whereas Eq. (|31|) for momentum 
wave packets fails. For a wave packet localized in momen- 
tum in Fig. lb), the momentum- wave-packet expression 
(|31fl works well and it is almost indistinguishable from the 
general DR lfTS|) . but the original position- wave-packet 
expression 1)27)) from Ref. ^| fails completely. The gen- 
eral DR works very well in both cases. It might seem that 
either the momentum or position versions could cover the 
whole range of Gaussian wave packets, because one might 
think that the intermediate case, i. e., a fairly symmet- 
ric wave packet, is localized enough in both position and 
momentum. That this is not so is provided by the fi- 
nal test in Fig. lc): both specialized expressions 1)30)) 
and 1)31)) give a significant error in comparison with ex- 
act fidelity, but the general DR l|18l) gives very accurate 
results, as expected because of its "fair" treatment of po- 
sition and momentum. To conclude, expression (|18|l . is 
accurate for the whole range of Gaussian wave packets, 
from position-like to symmetric to momentum-like, even 
in the presence of mixed dynamics. 
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B. Nonlocal states 

For nonlocal states, there is even less hope that the 
position-wave-packet expression for fidelity (|30|l from 
Ref. 13] would work. One might think that for a 
superposition of localized wave packets it is enough to 
simply add the terms H30|) for the fidelity amplitude. 
This is not the case which can be seen by considering 
a wave packet ip that is a superposition of two Gaussian 
wave packets ip\ and ip%, centered at phase space points 
(RijPjJ and (R 2 ,P 2 ). The resulting wave packet has 
a Wigner function that is not just a simple sum of the 
Wigner functions of the two Gaussian wave packets. The 
correct Wigner function has in addition an interference 
term localized in the vicinity of the phase-space point 
((Rj + R 2 )/2, (P 1 + P 2 )/2). We will demonstrate now 
the importance of this interference term and show that 
if it is taken into account, the general DR 118(1 will still 
give excellent results, even for nonlocal states. 

Being motivated by the quantum computation appli- 
cations, let us consider a superposition of computational 
states (i. e., position states in the abstract phase space), 
instead of Gaussian wave packets. Our initial state is a 
coherent superposition 

|V> = (|Ri) + |R 2 » , (35) 

with a Wigner distribution, 

p^ h (r, p) = i (2nh)- d {S(r - R x ) + 8(r - R 2 ) 
+ 26 [r - (Rx + R 2 ) /2] cos [(R x - R 2 ) • p/h]}. (36) 

If the interference term is neglected, we obtain a Wigner 
function of the incoherent superposition l|H8jl. 

Pw° h (r, P) = ~ (2^)~ d [S (r - R x ) + 6 (r - R 2 )] (37) 

Figure 2 compares two approximate ways to compute 
fidelity with the exact quantum result: both approxima- 
tions use the general DR i|18|) . but whereas one uses the 
correct full Wigner function (|36f) (black solid line), the 
other uses the incorrect Wigner function (|5Tjl , neglecting 
the interference term (purple dashed-dotted line). Again, 
the exact result is represented by solid dots. The param- 
eters used in Fig. 2 are n = 200, k = 0.7, e = 0.02, 
Qi = 0.47T, and 400 classical trajectories were used. The 
position of the other component state varies in the two 
parts. 

If the positions Ri and R 2 are largely separated, the 
oscillations in the interference term have a high fre- 
quency. Because nearby initial conditions follow similar 
trajectories and have similar actions, the phase factor 
in the DR (|18|) varies slowly. Therefore the fast oscilla- 
tions in the weight factor given by the interference term 
in the Wigner function can completely cancel out the 
contribution of the interference part to the DR integral. 
(Incidentally, this situation is in a way opposite to the 
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FIG. 2: Fidelity for a nonlocal state-coherent superposition 
of two position states-in a perturbed standard map with n = 
200, k = 0.7, e = 0.02. Comparison of the exact result (solid 
dots), general dephasing representation (black line) and the 
approximate DR neglecting coherence effects in the Wigner 
function (purple dashed-dotted line). The two position states 
are located at: a) Qi — OAn and Q2 = 1.2-7T, b) Qi — OAtt 
and Q2 = 0.427T. In both parts, 400 classical trajectories were 
used. 

usual semiclassical considerations where the weight is a 
slowly varying function and the phase factor is the fast 
oscillating factor.) Figure la) shows an example of sit- 
uation where this cancellation occurs: Q\ — 0.4-7T and 
Q2 = 1.2ir. Because the interference term is negligible, 
both approximations give the same and very accurate 
result. 

If the initial states are closer, as in Fig. 2b), where 
Qi = 0.47T and Q 2 = 0.42-7T, the interference term is 
important, and only the correct Wigner function l|36() 
agrees well with the exact result. This shows that for 
coherent nonlocal states, the general DR H18|) must be 
used instead of some approximate versions which neglect 
quantum coherence of the initial state. 

C. Mixed states 

Wigner function (|3T|) was wrong for the coherent state 
(|o3|l . but it does correctly describe a certain mixed state, 
namely the incoherent superposition of computational 
states |Ri) and |R 2 ), 

^ coh =i(|R 1 )(R x | + |R 2 )<R 2 |). (38) 

In Sec. IIIII it was shown that if the generalized defini- 
tion l|23|) of fidelity for mixed states is used, dephasing 
representation (|18|l remains valid, as long as the Wigner 
transform of the density operator (|26|l is used. For the 
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incoherent mixture with density operator (|38|l . Wigner 
distribution is precisely that given by Eq. 1(37(1 . Figure 
3 compares DR (118(1 with the Wigner function ((37(1 with 
the exact fidelity for the state ((38(1 . Fidelity computed 
by the DR is drawn with a black solid line, exact fidelity 
with solid dots. The parameters are the same as in Fig. 
2b), in particular Qi — OAir and Qi = 0.427T. Although 
now only 200 classical trajectories were used, the agree- 
ment is again excellent. 
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FIG. 3: Fidelity for a mixed state-incoherent superposition 
of two position states-in a perturbed standard map. All pa- 
rameters are the same as in Fig. 2b), except that only 200 
trajectories were used. Comparison of the exact result (solid 
dots) and general dephasing representation (black line). 

Last but not least we consider the completely ran- 
dom mixed state. It is an incoherent superposition of all 
computational states and in a finite-dimensional Hilbcrt 
space, its density operator is 

n 

p random = -V|Q i )(Qi| = -L 
n * — ' n 

i—l 

Figure 4 compares the random-state version 1(34(1 of DR 
(black solid line) with the exact result (solid dots). Pa- 
rameters in this calculation are n = 100, k = 2, e = 0.03 
and 1000 classical trajectories were used. Again, it is re- 
assuring that even in the case that the whole phase space 
is important, with just 1000 trajectories, dephasing rep- 
resentation still works so well-despite the fact that it was 
derived solely from semiclassical arguments and requires 
only classical information. 



VI. RELATION TO OTHER "WIGNER" 
METHODS 

It should be noted that the Wigner distribution has 
been used in various other approximate methods, es- 
pecially in chemical physics. For instance, it was used 
to compute photodissociation cross-sections (53, |53 , to 
treat inelastic scattering [6(j , or to compute thermal cor- 
relation functions using the linearized semiclassical IVR 
method [44l fell 62]. In all these applications, there was 
just one Hamiltonian, but the two states (or more gen- 
erally, density or other operators) were different. The 




10 20 30 40 50 
t 



FIG. 4: Fidelity for a random mixed state-incoherent super- 
position of all basis states-in a perturbed standard map with 
n = 100, k — 2, e — 0.03. Comparison of the exact (solid dots) 
result and general dephasing representation (black line). 1000 
classical trajectories were used. 

quantity of interest was a general correlation function of 
the type 

C A B{t) = ti(AU^BU) (39) 

where A and B are general operators and U = Te~* J HdT 
is the time evolution operator. Usin g va rious approxima- 
tions, all authors 

El, EE E3, E3,l 

1621 obtain the same 
final result, expressed as an overlap of two Wigner dis- 
tributions, one at time 0, the other evolved classically to 
time t, 

^Wigner (fi = ( 2ir H) d [ drQ J dpoAw (r ,p )B w (r t ,p t ) 

(40) 

Here Aw and Aw are the Wigner transforms ((261) of 
operators A and B. 

Because there is only one Hamiltonian, there is no de- 
phasing factor e lAS / h : as in the DR. In fact we could 
apply one of these older approaches to the second gener- 
alized definition 1(241) of fidelity for mixed states because 
that definition is in the form of Eq. ((391) with A = B = p 
and the time evolution operator U = e +lH t / h e - lH l / h _ 
Then we would obtain a very different result from the 
DR, 

M Wigncr {t) = {2nh)d J dj . Q J dpQf)w (ro5 po) pw 

(41) 

Although appearing as elegant as the dephasing represen- 
tation, there is a problem with this expression. First, it 
will be much more sensitive to numerical errors. We can 
see that already by considering zero perturbation. Cor- 
rectly, for each initial condition ro , po , we should have 
ro = r f and p t = po. In systems with nonlinear dy- 
namics, particularly chaotic systems, numerical errors in 
forward and backward propagation will yield exponen- 
tially growing errors. If the initial state is a localized 
wave packet, expression ((4111 would give a numerically 
decaying overlap even for zero perturbations when exact 
fidelity is constant M(t) = 1. Indeed, numerical test not 
presented here showed that instead of staying at unity, 
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M Wlgncr quickly decays to a plateau and remains there 
for some time, and finally decays exponentially again. 
(This is the same behavior as observed in literature for 
physical perturbations 0, 0] . 

Even if numerical errors did not exist, equation i|41|) 
would have problems. It can describe some decay due 
to dephasing, but only that in the fast oscillating parts 
of the initial state. For simple Gaussian wave packets, 
the fidelity decay in Eq. I|41|) is completely due to the 
decay of classical overlaps, i. e., classical fidelity. To con- 
clude, the "Wigner" form (|41|l is apparently not as good 
as the dephasing representation, but it does deserve fur- 
ther study, especially because it might shed further light 
on the question of importance of various contributions 
to fidelity. Preliminary studies show that M Wlgner cor- 
rectly describes exact fidelity in both chaotic and quasi- 
integrable systems for large perturbations (i. e., in Lya- 
punov and algebraic regimes, respectively), when dephas- 
ing is not important |15|. It gives wrong results in both 
chaotic and quasi-integrable systems for small perturba- 
tions (in the FGR and Gaussian regimes), when dephas- 
ing is important [l5| . 

VII. CONCLUSION 

This paper has presented a derivation of a general semi- 
classical expression for fidelity of pure and mixed states. 
This dephasing representation expresses fidelity as an in- 
terference integral, with weight of each term given by the 
Wigner function and the phase by the integrated pertur- 
bation along an unperturbed trajectory. In particular, 
no analog of the Van Vleck determinant is needed. As 
the original specialized expression llSU|l from Ref. [l3| . 
dephasing representation avoids searching for the expo- 
nentially growing number of terms in the standard semi- 
classical expressions 3] • It also avoids the ubiquitous di- 



vergences in Van Vleck determinants present in the usual 
semiclassical expressions. 

The advantage of dephasing representation lies in that 
it does not require the original state to be localized. Its 
form suggests that it should be applicable to general pure 
and mixed states. This claim was supported by the fol- 
lowing numerical evidence: First, it was shown, on the 
example of Gaussian wave packets, that position and mo- 
mentum must be treated symmetrically. This was the 
flaw of the expression from Ref. isj and is apparently 
corrected in the DR. Second, on the example of coher- 
ent superpositions of states, it was shown that oscillatory 
patterns in the Wigner function are important: therefore 
classical phase space distribution, resulting from inco- 
herent superposition of component Wigner distributions 
(for states for which these are the same as classical dis- 
tributions). This may shed some further light on the 
controversial issue of importance of sub-Plank structures 
on decoherence |22|, |63j . Finally, it was shown that DR is 
also accurate for mixed states: incoherent superpositions 
and completely random states. All tests were performed 
on a system with mixed phase space: with both inte- 
grable and chaotic regions. 

While the numerical tests were quite successful, a fur- 
ther study is needed to determine precisely all situations 
where the dephasing representation breaks down. The 
analysis provided in Sec. [n]of this paper should simplify 
that task. Also, a more rigorous formulation of the pre- 
cise conditions of validity of the dephasing representation 
is needed. 
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